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Abstract 

The cross section of the process e'’'e“ —> 7r'’“7r“ has been measured using about 
114000 events collected by the CMD-2 detector at the VEPP-2M e+e” collider 
in the center-of-mass energy range from 0.61 to 0.96 GeV. Results of the pion 
form factor determination with a 0.6% systematic uncertainty are presented. The 
following values of the p- and (u-meson parameters were found: Mp = (776.09 ±0.81) 
MeV, Tp = (144.46 ± 1.55) MeV, r(p ^ e+e") = (6.86 ± 0.12) keV, Br{u; 

= (1.33 ± 0.25)%. Implications for the hadronic contribution to the muon 
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anomalous magnetic moment are discussed. 


1 Introduction 


The cross-section of the process e’''e —7r’''7r is usually expressed in terms of the pion electromag¬ 
netic form factor Ft^{s): 




ttq; 




( 1 ) 


where s is the center-of-mass (c.m.) energy squared, is the pion mass and = y 1 — is 

the pion velocity in the c.m. frame. 


The pion form factor measurement is crucial for a number of physics problems. Detailed experimental 
data in the time-like region allow a determination of the parameters of the p(770) meson and its 
radial excitations. Extrapolation of the energy dependence of the pion form factor to the point s = 0 
provides a value of the pion electromagnetic radius. 

In the energy range below 0(1020) the process e’''e“ —tt+tt” gives the dominant contribution to the 
quantity R{s) dehned as 

R = a{e~^e~ hadrons)/a{e'^e~ jjT). 

R{s) is an important measurable quantity widely used for various QCD tests as well as calculations 
of the dispersion integrals. For such applications, at high energies R{s) is usually calculated within 
the perturbative QCD frame, while for the low energy range the direct measurement of the e’''e“ — 
hadrons cross-section is necessary. 


Particularly, knowledge of R{s) with high accuracy is required for the evaluation of the hadronic 
contribution to the anomalous magnetic moment of the muon (g-2)/2 (see [1] and references 
therein). The ultimate goal of the experiment E821 [2] running in Brookhaven National Laboratory is 
to measure the muon anomalous magnetic moment with a relative precision of 0.35 ppm. Within the 
Standard Model (SM) the uncertainty of the theoretical value of the leading order a^ is dominated 
by the uncertainty of the hadronic contribution calculated via the dispersion integral 
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where -R(s) is the QED kernel and (t(s) is the cross section of hadrons. The precision of the 

ahad calculation depends on the approach used and varies from 1.34 ppm based on data only [3] 
to 0.53 ppm if in addition r-lepton decay data as well as perturbative QCD and QCD sum rules are 
extensively used [4]. The major contribution to its uncertainty comes from the systematic error of 
the R(s) measurement at low energies (s < 2 GeV^), which, in turn, is dominated by the systematic 
error of the measured cross section e+e“ —tt+tt”. 


Assuming conservation of the vector current (CVC) and isospin symmetry, the spectral function of 
the T~ —> decay can be related to the isovector part of the pion form factor [5]. The detailed 
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measurement of the spectral functions was provided by ALEPH [6], OPAL [7] and CLEO-II [8]. The 
comparison of the pion form factor measured at colliders with the spectral function of the 

T~ —> decay provides a test of CVC. If CVC holds with high accuracy, r-lepton decay data 

can be also used to improve the accuracy of the calculations mentioned above [9,10]. 

E821 has recently published the result of its measurement of with an accuracy of 1.3 ppm [11]. 
The measured value of is 2.6 standard deviations higher than the SM prediction of [4]. Q This 
observation makes new high precision measurements of the —>■ hadrons cross section and 

particularly of the pion form factor extremely important. 

Since early 70s the VEPP-2M collider has been running in the Budker Institute of Nuclear 

Physics in the c.m. energy range 360-1400 MeV. The most precise pion form factor data were obtained 
in late 70s ~ early 80s by CMD and OLYA detectors [13]. Their accuracy was limited by systematic 
errors of the experiments, varying from 2% to 15% over the VEPP-2M energy range. 

The CMD-2 detector installed in 1991 is a general purpose detector consisting of the drift chamber, 
the proportional Z-chamber, the barrel Csl calorimeter, the endcap BGO calorimeter installed in 
1996, and the muon range system. The drift chamber, Z-chamber and the endcap calorimeters are 
placed inside a thin superconducting solenoid with a held of 1 T. More detail on the detector can be 
found elsewhere [14]. 

Though the collider luminosity has not considerably increased since the previous work in the 1980’s, 
the present detector design is signihcantly improved, resulting in a smaller systematic error of the 
pion form factor in this measurement. Particularly, the following advantages of CMD-2 should be 
mentioned compared to previous detectors: simultaneous measurement of the particle momentum 
and energy deposition which simplihed particle identihcation and helped to reduce background; high 
precision determination of the hducial volume with the help of the Z-chamber; a thinner beam pipe 
which reduced the nuclear interactions of pions. 

During 1994-1995 a detailed scan of the c.m. energy range 610-960 MeV was performed. Since the 
pion form factor changes relatively fast in this energy range dominated by the p-meson, the systematic 
error due to an uncertainty in the energy measurement can be signihcant. To reduce this contribution 
to a negligible level, the beam energy was measured with the help of the resonance depolarization 
technique with an accuracy of 140 keV for almost all energy points [15]. This paper presents the 
hnal analysis of the data taken in 1994-1995 with the integrated luminosity of about 310 nb“^ and a 
systematic uncertainty of the cross section of 0.6%. Preliminary results with a systematic uncertainty 
of 1.4% were published in [16], based on the same data sample. 


2 Data Analysis 


The data were collected at 43 points with c.m. energy ranging from 610 MeV to 960 MeV in 10 MeV 
energy steps, except for the narrow energy range near the a;-meson, where the energy steps were 2-h6 
MeV. 

^ Recent progress in estimating the light-by-light scattering contribution to [12] implies that the dif¬ 
ference between experiment and theory reduces to about 1.5 standard deviations. 
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Fig. 1. Spatial distribution of the vertex. The left plot shows the distribution of the distance from the vertex 
to the beam axis {p), the right plot presents the distribution of the distance between the vertex and the 
center of the detector along the beam axis {z). The open histogram corresponds to all collinear events, the 
filled one shows the subset of background events. 

From more than 4 ■ 10^ triggers recorded, about 3 ■ 10^ events were selected as collinear, with a 
signature of two particles of opposite charge and nearly back-to-back momenta originating from the 
interaction point. The following selection criteria were used: 

(1) Two tracks of opposite charge originating from the interaction region are reconstructed in the 
drift chamber. 

(2) The distance from the vertex to the beam axis, p, is less than 0.3 cm and the z-coordinate of 

the vertex (along the beam axis) is within —15 < < 15 cm. 

(3) The average momentum of the two particles (pi +p2)/2 is between 200 and 600 MeV/c. 

(4) The difference between the azimuthal angles (in the plane perpendicular to the beam axis) of 
two particles lA^?! = \ti — \(pi — (p 2 \ \ < 0.15. 

(5) The difference between the polar angles (the angle between the momentum and the beam axis) 
of two particles |A0| = |0i — (tt — 02)| < 0.25. 

(6) The average polar angle of two particles Qavr = [0i + (7J'—02)]/2 is within 1.1 < 0^^^ < (vr —1.1). 
This criterion determines the hducial volume. 

The selected sample of collinear events contains 7r’''7r“, p~ 

events (below referred to as beam originating) as well as the background of cosmic particles which 
pass near the interaction region and are misidentihed as collinear events. The number of cosmic 
background events Ncosmic was determined by the analysis of the spatial distribution of the vertex. 
Both distributions of the longitudinal coordinate {z) and the distance from the beam axis (p) are 
peaked around zero for the beam originating events, but are very broad, almost flat for the cosmic 
background events. Typical p- and ; 2 -distributions are shown in Fig. 1. 

The momenta of e, p, tt from the beam originating events are rather close (the difference is comparable 
to the momentum resolution of the drift chamber), so the overlap of the momentum distributions is 
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Fig. 2. Energy deposition of collinear events for the beam energy of 400 MeV. 
large. Therefore, separation of e, p and tt by their momenta is impossible at > 600 MeV. 

On the contrary, the energy deposition of the particles in the calorimeter is quite different for e, p 
and TT. The typical energy deposition of two particles {E^ vs E~) for experimental events selected 
at the beam energy of 400 MeV is shown in Fig. 2. The high deposition spot corresponds to e+e“ 
pairs, where both particles leave almost all their energy in the calorimeter. The low deposition spot 
represents pairs, cosmic muons and those 7r’''7r“ pairs, in which both particles interact as 

minimum-ionizing. The long tails correspond to 7r’''7r“ pairs in which one or both particles undergo 
nuclear interactions inside the calorimeter. 

Therefore, the energy deposition of the particles was used for the separation of the beam originating 
events. Since the overlap of the distributions for e+e“ and 7r’''7r“ pairs is small, this approach gives 
stable results with a small systematic error. However, and tt+tt" pairs cannot be separated 

well by their energy deposition. To avoid this problem, the number of pairs was derived from 

the number of e+e” pairs according to QED, taking into account radiative corrections and detection 
efficiencies. Since in this energy range the number of pairs is small compared to that of 

the systematic error caused by the corresponding calculation is negligible (less than 0.03%). 

The separation was based on the minimization of the following unbinned likelihood function: 

i = - E In (E ■ UE*, E-)] + ]V„ (3) 

events \ cl /a 

where a is the event type (a = ee, p/i, tttt, cosmic), Na is the number of events of the type a and 
fa{E~^,E~) is the probability density function (p.d.f.) for a type a event to have energy depositions 
E^ and E~. It was assumed that E~^ and E~ are not correlated for events of the same type, so the 
p.d.f. can be factorized as 


UE+,E-) = f+{E+)-f-{E-), 

where fa{E) is the p.d.f. for e^, and cosmic muons to have energy deposition equal to E. This 
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assumption is not entirely correct since there are small correlations between E~^ and E~ because of 
the dependence of the calorimeter thickness on the polar angle as well as the c.m. energy shift due 
to the initial state radiation. The hrst effect was corrected for while the second one was studied with 
the help of simulation and was shown to be negligible (below 0.1%). 

For pairs and cosmic events the energy deposition does not depend on the particle charge, 

while the energy depositions for and 7i~ are different. Therefore = f~ for a = ee, /r/i and 
cosmic, but in case of pions the corresponding p.d.f.’s are different. 

As was mentioned before, the ratio N^^/Nee was hxed during minimization according to the QED 
calculation 

(1 + 

Nee ^ee (1 + ^ee)^ee 

where a is the Born cross section, 5 is the radiative correction and e is the detection efficiency which 
includes acceptance as well as reconstruction and trigger efficiencies. The number of cosmic events 
Ncosmic determined separately was hxed during the minimization while its huctuation was added in 
quadrature to the huctuation of 

To obtain the specihc form of p.d.f.’s, the energy deposition of e, fi and tt in the CMD-2 calorimeter 
was studied. For electrons (positrons) and cosmic muons it can be obtained with the help of the data 
itself. Particles of positive charge with large enough energy deposition are almost 100% positrons and 
they were used to tag electrons. Events with a large value of p are mostly cosmic muons. Such tagged 
particles were used to determine the energy deposition of electrons and cosmic muons respectively. 

The simulation was used to obtain the energy deposition of muons from e’''e“ — pNp,~. In the 
energy range under study these muons interact purely as minimum-ionizing particles which are well 
described by the simulation. 


On the contrary, the simulation of the interaction of low energy pions inside the calorimeter is 
not reliable enough. In addition, there is no good method to tag pions from e+e“ —tt+tt" events. 
Therefore, the p.d.f.’s for and 7r“ were obtained from the analysis of the energy deposition of pions 
coming from the 0(1020) —decay. Erom a large data sample collected by CMD-2 around 
the 0-meson peak, about 10^ 0(1020) —> Stt events were selected with practically no background. 
Pions found in these events cover the whole interesting range of momenta and angles. The energy 
deposition of these pions was analyzed and the parameterization for tt’*' and 7r“ p.d.f.’s was derived. 


Finally, to simplify the final error calculation, the likelihood function (3) was rewritten to have the 
following global fit parameters: 


[Nee + N^^), 




Nee + N, 




instead of Nee and N^^t, (with N^^/Nee and Ncosmic fixed). The likelihood function has some other 
fit parameters characterizing the p.d.f.’s for different types of particles, such as the mean energy, 
energy resolution, the asymmetry of some distributions etc. More detail about the energy deposition 
of different types of collinear events and the p.d.f. parameterization can be found in [16]. 


After the separation procedure the following number of events was obtained for three described above 
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Fig. 3. Difference between the initial and the reconstructed values of the form factor for the simulated data 
set. 

classes: N^e + = 180038, = 113824 and Ncosmic = 17390. 

Special studies were performed to estimate the systematic error of the separation procedure. The 
dominant effect was produced by the small nonuniformity of the calorimeter calibration. Due to the 
forward-backward asymmetry of the e+e“ —> e’''e“ cross section, the calorimeter calibration error 
leads to a small difference between and e~ energy depositions. The corresponding error was found 
to be less than 0.2%. Several different functional forms were used to parameterize p.d.f.’s of e’s and 
tt’s and the hnal cross section was stable within 0.1% for different selections. The existing variation 
of the calorimeter response between calibrations leads to small variations of the energy resolution 
which could also influence the results. The estimated contribution of this effect is below 0.1%. As a 
hnal test, the large amount of > e’''e“( 7 ), events was generated in a proper 

proportion with the help of full detector simulation at several energy points covering the whole energy 
range. After that the simulated events were subject to the same separation procedure as for the data. 
Results are shown in Fig. 3. The reconstructed value is always consistent with the input one and the 
average difference is below 0.1% or consistent with zero. Thus, the systematic error because of the 
event separation taking into account the above effects is estimated to be 0.2%. 


3 Form factor calculation 


The pion form factor was calculated as: 

1^ |2 _ -^TTn ^ ^ee ‘ (1 T Ae)^ee T ^' (7 T ^^ 

A^ee + O’TTTr ■ (1 + <^7r7r)(l + Ajv)(l -|- 

where the ratio + was obtained from the minimization of (3), a are the corresponding 

Born cross-sections, 5 are the radiative corrections, e are the detection efficiencies, Aj) and Ajv are 
the corrections for the pion losses caused by decays in flight and nuclear interactions respectively, 
and Astt is the correction for misidentihcation of cn —events as the e+e“ —tt+tt". In the 
case of the e+e“ —tt+tt" process, corresponds to pointlike pions. 

Radiative corrections were calculated according to Refs. [17,18] in which the accuracy of the obtained 
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formulae was 0.2%. The improved precision of the calculations compared to previous works [19] 
comes from taking into account the radiation of the additional photon in a narrow cone along the 
direction of electrons and positrons. The radiative corrections were calculated by the Monte Carlo 
integration of the differential cross sections imposing all selection criteria. We estimate the uncertainty 
of the form factor because of the radiative corrections to be 0.4% dominated by the accuracy of the 
ratio (1 + 5ee)/(l + 

The radiative corrections for the process e’''e“ —>■ 7r’''7r“ include the effects of both initial (ISR) and 
hnal state radiation (FSR) and do not include the vacuum polarization terms (both leptonic and 
hadronic) since the latter are considered to be an intrinsic part of the hadronic cross section and 
corresponding form factor (1). However, for various applications based on dispersion relations and 
involving the total cross section of e’''e“ —hadrons, the radiation by hnal pions is no longer a radiative 
correction, so that with a photon radiated by one of the hnal pions should be considered as 

one of the possible hadronic hnal states contributing to the total cross section. Therefore, the bare 
cross section e+e“ —>• 7r+7r“(7), below referred to as was also calculated as 

y.,,.) = |F,(i>)|=. |1 - n(i,)|^. (l + . (5) 

The factor |1 — n(s)p with a polarization operator n(s) excludes the ehect of leptonic and hadronic 
vacuum polarization, so that one obtains the bare cross section required for various applications. 

A correction for the hnal state radiation A(s) was calculated based on [18], where only the ehects of 
FSR, integrated over the whole allowed kinematic region, were taken into account. The result 


A(S) = i±a{4Li2(i^) +2Li2(-ij|i) - |31n(j^) + 21n;3.1L^,} 


■1-A 


l-/3n 


■ 1+A 


-SNlTf)-4111/3.,+ i 4(l+®=-2 




l + Z^TT 
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' l + Z^TT 
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where = In Li 2 (^) = - Jo f ln(l - x), coincides with [21,20]. 


( 6 ) 


Two diherent trigger settings were used during data taking. For energies below 810 MeV only a 
charged trigger was used in which the positive decision was based on the information from the track¬ 
ing system only. The efficiency of the charged trigger was measured to be higher than 99% and, 
more important, equal for diherent types of collinear events in the energy range under consideration. 
Therefore, the uncertainties related to the trigger efficiency cancel in (4). Above 810 MeV the addi¬ 
tional requirement that the energy deposition in the calorimeter be at least 20 MeV was applied for 
the trigger. The efficiency of this requirement was measured to be 99.5% for e’''e“ —tt+tt", 99.2% for 
/i+yU" and 100% for e+e“ —>• events, and was included in the corresponding detection 
efficiency. These corrections give a negligible contribution to the systematic uncertainty. 


The reconstruction efficiency was measured using the experimental data themselves. It was found to 
be within the 98-100% range at all energies and, within the statistical accuracy, nearly the same for 
all types of collinear events. Therefore, it cancels in (4). The systematic error of the cancelation was 


^ We corrected a misprint in formula (2.5) of Ref. [18], also noted by the authors of Ref. [20]. 

^ Note that the definition of the cross section above is different from that in [16] where the FSR cross 
section was not taken into account and only the correction for the leptonic vacuum polarization was applied. 
Therefore, their comparison is meaningless. 













estimated to be better than 0.2%. It’s worth noting that such a cancelation allows a determination 
of the form factor with better precision than that of the luminosity. 

The hducial volume (detection solid angle) is determined by selecting events with the average polar 
angle Oavr = (©i + tt — 02)/2 in the range between ©mm and (tt — ©min) with ©min = 1.1. The value 
of Oavr is determined by the CMD-2 Z-chamber [22] which has a spatial resolution along the beam 
axis better than 0.6 mm. That corresponds to a systematic error in the form factor of about 0.2%. 
To test this estimate, the pion form factor was also determined at ©min = 1-0 radian. The difference 
between the form factor values determined at two values of ©min averaged over the c.m. energy range 
was found to be (0.1 ± 0.3)%, consistent with zero. 

The correction Ad for the pion losses caused by decays in flight was calculated with the help of 
simulation. Its value, varying from 0.2% at 600 MeV to 0.03% around the 0(lO2O)-meson, turned 
out to be small because of the small size of the decay volume and the fact that the maximum decay 
angle is small enough and of the same magnitude as the angular resolution of the drift chamber. This 
effect gives a negligible contribution to the systematic uncertainty. 

The correction Ajy for the pion losses caused by nuclear interactions inside the wall of the beam pipe 
or the drift chamber was calculated with the help of FLUKA-based simulation [23]. The value of 
the correction is slowly changing from 1.7% at 600 MeV to 0.8% at the (;/)(1020)-meson energy. The 
systematic error of 0.2% for the A^r calculation was estimated from the uncertainty of the nuclear 
cross sections used in FLUKA [24]. 

There is also a small correction for the losses of > e’''e“ events, caused by interactions of 

electrons with material of the beam pipe and the drift chamber. It was taken into account by reducing 
A AT by 0.15% according to the simulation. 

Use of the resonance depolarization for the absolute beam energy calibration reduced the systematic 
uncertainty of the form factor due to the beam energy measurement to 0.1%. 

In the narrow energy range around the ci;(782)-meson there is a small background of e’''e“ — 
events, misidentihed as 7r’''7r“. The corresponding correction As^ was calculated from the 

simulation. It reaches its maximum value of about 1% at the cc;(782)-meson energy and drops fast to 
nearly zero at the energies outside the cu-meson. 

The measured values of the pion form factor as well as those of the bare 7r’''7r“(7) cross- 

section obtained from (5) are shown in Table 1. Only statistical errors are given. The main sources 
of the systematic error are listed in Table 2. The overall systematic error obtained by summing 
individual contributions in quadrature is about 0.6%. The values of the pion form factor and the 
bare cross section in Table 1 supersede our preliminary results presented in Table 3 of Ref. [16]. 


4 Fit to data 


The parameterization of the pion form factor in the energy range under study should include contri¬ 
butions from the p(770), ci;(782) and p(1450) resonances. It was assumed that the only mechanism 
for the LO 7r~^7r~ decay is p — a; mixing. Following [25], we represent the wave function of the 
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Table 1 

The measured values of the pion form factor and bare cross section e^e“ ^ 7r'’“7r“(7). Only statistical 


errors are shown. The systematic error is estimated to be 0.6%. 


2E (MeV) 

to 

/ N, nb 

7r7r(7) ’ 

2E (MeV) 

\FfI 


610.50 

8.00 ± 1.13 

327.8 ± 46.2 

784.24 

34.30 ± 0.97 

935.9 ± 26.5 

620.50 

9.71 ±0.72 

390.6 ± 29.1 

786.04 

30.01 ± 1.07 

809.5 ± 28.9 

630.50 

10.82 ± 0.73 

426.9 ± 28.6 

790.10 

31.96 ± 1.11 

859.3 ± 30.0 

640.51 

11.02 ±0.72 

426.2 ± 27.9 

794.14 

29.81 ±0.85 

799.6 ± 22.7 

650.49 

12.66 ±0.83 

480.3 ±31.6 

800.02 

29.40 ± 0.68 

782.7 ± 18.1 

660.50 

13.40 ± 0.75 

498.8 ± 27.8 

810.14 

25.60 ±0.55 

669.7 ± 14.4 

670.50 

15.36 ±0.82 

560.5 ± 30.0 

820.02 

24.09 ± 0.78 

619.0 ±20.0 

680.59 

19.49 ±0.90 

697.2 ± 32.3 

829.97 

20.52 ± 0.73 

517.5 ± 18.5 

690.43 

20.49 ± 0.73 

718.6 ±25.5 

839.10 

16.75 ±0.70 

415.4 ± 17.3 

700.52 

23.81 ±0.62 

818.1 ±21.2 

849.24 

14.16 ±0.70 

344.8 ± 17.0 

710.47 

27.52 ± 0.99 

926.1 ±33.4 

859.60 

14.47 ± 0.67 

345.7 ± 15.9 

720.25 

31.80 ±0.82 

1048.1 ±27.0 

869.50 

11.10 ±0.46 

260.7 ± 10.7 

730.24 

34.34 ± 1.18 

1107.2 ±38.1 

879.84 

10.13 ±0.79 

233.6 ± 18.1 

740.20 

37.82 ± 1.18 

1191.9 ±37.1 

889.72 

8.44 ± 0.34 

191.1 ± 7.7 

750.28 

42.33 ± 1.18 

1302.9 ± 36.2 

900.04 

7.74 ± 0.29 

172.1 ±6.4 

760.18 

43.44 ± 1.21 

1306.9 ± 36.5 

910.02 

6.82 ±0.31 

149.1 ±6.8 

764.17 

43.23 ± 1.03 

1289.6 ± 30.6 

919.56 

6.03 ±0.30 

129.7 ±6.5 

770.11 

44.23 ± 1.14 

1304.6 ± 33.5 

930.11 

5.70 ± 0.36 

120.3 ±7.7 

774.38 

43.19 ± 1.15 

1265.7 ±33.7 

942.19 

5.16 ±0.25 

106.8 ±5.1 

778.17 

45.41 ± 1.29 

1321.6 ±37.7 

951.84 

4.56 ±0.23 

92.9 ±4.8 

780.17 

42.90 ± 1.21 

1234.9 ± 34.8 

961.52 

4.29 ±0.23 

86.2 ±4.6 

782.23 

36.59 ±0.66 

1026.6 ± 18.5 





ci;(782)-meson as 


— |<^o) + ^Ipo) 5 


where |a;o) and |po) are the pure isoscalar and isovector respectively, and e is the p — ui mixing 
parameter. Then, in the energy region close to the p(770)- and a;(782)-meson masses, the form factor 
can be written as 


F^{s) 


F F 

p +e ^ 
s-Mj s-Ml 


Fp 

-Ml 


+ e 


-M2 



1 + ^ 


FJM^ - ■ 

- AQ) J ’ 


(7) 


where we keep only the terms linear in e. The quantities and Mp are complex and contain the 
corresponding widths. 
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Table 2 

Main sources of the systematic errors 


Source 

Contribution 

Event separation 

0.2% 

Radiative corrections 

0.4% 

Detection efficiency 

0.2% 

Fiducial volume 

0.2% 

Correction for pion losses 

0.2% 

Beam energy determination 

0.1% 

Total 

0.6% 


Combining together the contributions from the p(770)- and p(1450)-mesons, and including that from 
the ci;(782)-meson as in (7), we write the pion form factor: 


F^{s) = 




1 + 5 


M; 


BW. 



+ /5-BW«S45o)(^) 


1 + /3 


( 8 ) 


where parameters S and (3 describe the contributions of the ci;(782)- and p(1450)-mesons relative 
to the dominant one of the p(770)-meson. For the p(770) and p(1450) the GS parameterization is 
used [26]: 


BW 




M2(i + d.r,/M,) 
M2-s + /(s)-*M,r,(s)’ 


where 


(9) 


f{s) = r 


iVJ.p 




pI{s) (h{s) - + (Mj - s)pI{MI) ^ 


ds 


s=Ml 


( 10 ) 


h{s) 


yi + 2p^(s) 

- 1 =- In- 

TT ^/s 2mn 


and d is chosen to satisfy BW^(^^)(0) = 1: 


d = 


3 ml M, + 2p^{Ml) 

7rp2(M2) 2m^ 


+ 


Mr 


27rpM] 


mlMp 

7rp3(M2)’ 


( 11 ) 


( 12 ) 


where Mp, Fp are the p(770)-meson mass and width, and Pn{s) is the pion momentum at the squared 
c.m. energy s. For the energy dependence of the p(770) width, the P-wave phase space is taken: 


U(i>) 


r. 


pA^) 

3 

[M^l 

±\j.p 

U(MJ)J 


s 


(13) 


For the a;(782)-meson contribution the simple Breit-Wigner parameterization with a constant width 
is used. 
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To obtain the p(770)-meson leptonic width r(p — e^e ), the well known VDM relations were used 
|27]: 


V- 


47ra^ 

3Mf. 


9Vy 


and r 


V- 


ffk. PUM^v) 

GTT My 


Assuming 


Qp’yQp'irn 


Mlii + d- r,/M, 


(1 + p) 


and r ^_^ 7 r“t" 7 r~ — T 


pi 


the following result was obtained: 


2aYn (m^ (1 + d ■ Tp/Mpf 
QMpTp (l + /3)2 


Similarly, assuming (14) and 


d-MMBW,%o)(MS)l 


QLwyQu. 


(1 + / 9 ) 

the following expression is obtained for the branching ratio of the u 7i~^7r~ decay: 


Br\(jJ TT^TT ) = 


QM r + -T 

02 —>e+e 02 




{I+Pf 


(14) 


(15) 


(16) 


(17) 


(18) 


It should be mentioned that our parameterization of the a;(782)-meson contribution is slightly dif¬ 
ferent from that in [6,13]: (1 -|- <5 ■ s/M"^ ■ BWo 2 (<s)) instead of (1 -|- <5 ■ BWo 2 (<s))/(l -|- d). Fits to either 
parameterization give the same result. 

To estimate the model dependence of the obtained parameters, another model of the form factor 
parameterization was considered - the Hidden Local Symmetry (HLS) model [27,28]. In this model 
the p(770)-meson appears as a dynamical gauge boson of a hidden local symmetry in the non-linear 
chiral Lagrangian. The p(1450) contribution is not taken into account, replaced by the non-resonant 
coupling 77r’''7r“. This model introduces a real parameter a related to this non-resonant coupling. 
The original parameterization of the pion form factor was modihed in a similar way to take into 
account the ci;(782)-meson contribution. 


5 Results and discussion 


5.1 p(770)- and u{782)-meson parameters 


There are several parameters of both models whose values have to be taken from other measurements 
[29,30]. Their values were allowed to fluctuate within the stated experimental errors. The following 
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Fig. 4. Fit of the CMD-2 pion form factor data to the GS model. 

values of the a;(782)-meson parameters were taken from the CMD-2 experiment [30]: = (782.71 ± 

0.08) MeV, = (8.68 ± 0.24) MeV, F^ee = (0.595 ± 0.017) keV. Parameters of the p(1450) were 
taken from [29]: Mp(i 45 o) = (1465 ± 25) MeV and rp(i 45 o) = (310 ± 60) MeV. 

Results of the £t in the GS model are shown in Fig. 4. Parameters of GS and HLS models obtained 
from the fit are in good agreement with each other and are shown in Table 3. The hrst error is 
statistical, the second one is systematic. Two effects were taken into account in the estimation of 
the latter: the systematic uncertainty of the form factor measurement of 0.6% (Table 2) and the 
contribution of the p(1700) resonance missing in the adopted GS parameterization. The parameter 
f3 effectively describes the overall contribution of the p(1450) and p(1700) and its value is strongly 
correlated with the p(1450) width. To extract the p(1450) contribution, a fit of the data should be 
performed in a broader energy range. As a result, the value of (3 is strongly model dependent and is 
not well defined in the energy range under study. For this reason, its systematic error is not shown. 

In Table 4 the final results obtained in the GS model are compared to the world average values 
[29]. The values of the p(770) mass and width shown there are based on the previous measurements 
at VEPP-2M [13], ALEPH [6] and CLEO [8]. While the mass of the p(770) obtained in this work 
is in very good agreement with the previous measurements, our value of the width is 2.7 standard 
deviations lower. This difference can partly be explained by the difference of our parameterization 
compared to the previous works. Particularly, the value of Pp is correlated with the value of arg 5. In 
previous papers the parameter 5 was assumed to be real. Fixing arg 5 = 0 increases our value of Fp by 
2 MeV. Such model uncertainties as the effect of the complex phase of 5 and the difference between 
the results of GS and HLS fits were not included into the systematic error of the fit parameters. 

The leptonic width of the p(770) is in good agreement with the result of [13] quoted by [29]. Our 
value of the branching ratio uj is 1.6 standard deviations lower than the world average 

(2.1 ± 0.4)% based on the two most precise measurements from e’''e“ experiments [13,31]. Again, 
a different parameterization of the form factor was used in previous works. Also note that the 
parameters of the ci;(782)-meson, such as mass, width and the leptonic width, have changed, that 
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Table 3 

Results from fits to |F 7 r(s)p for the GS and HLS models 


Parameter 


GS model 

HLS model 

Mp, MeV 


776.09 ± 0.64 ± 0.50 

775.23 ±0.61 ±0.50 

Tp, MeV 


144.46 ± 1.33 ± 0.80 

143.88 ± 1.44 ± 0.80 

T{p e+e“), 

keV 

6.86 ±0.11 ±0.05 

6.84 ±0.12 ±0.05 

Br{uj 7:~^7r~ 

-),% 

1.33 ±0.24 ±0.05 

1.32 ±0.24 ±0.05 

l<^l 


(1.57 ±0.15 ±0.05) • 10-3 

(1.57 ±0.15 ±0.05) • 10-3 

arg(5 


12.6° ±3.7° ±0.2° 

13.0° ±3.7° ±0.2° 

[3 (GS) 


-0.0695 ± 0.0053 

— 

a (HLS) 


— 

2.336 ± 0.016 ± 0.007 

3d tv 


0.92 

0.94 

Table 4 

Comparison of the fit to F^(s) p in the GS model to the world average values 


Parameter 


GS model 

World Average 

Up, MeV 


776.09 ± 0.64 ± 0.50 

775.7 ±0.7 

Pp, MeV 


144.46 ± 1.33 ± 0.80 

150.4 ± 1.6 

r(p^ e+e-), 

keV 

6.86 ±0.11 ±0.05 

6.77 ±0.32 

Br(u) TT~^TT~ 

-),% 

1.33 ±0.24 ±0.05 

2.1 ±0.4 


also affects the extracted value of the branching ratio. Our fit to the old data [13] gives the value 
Br{uj 7i^n~) = (2.00 ± 0.34)%, still 1.6 standard deviations above our present result. 


5.2 Hadronic contribution to the muon anomalous magnetic moment 


Let us estimate the implication of our results for the contribution from the annihilation into 
two pions, which dominates the hadronic contribution to (g-2)/2. To this end we compare its value 
in the energy range studied in this work and calculated from CMD-2 data only to that based on 
the previous e+e“ measurements [13,31]. Table 5 presents results of the calculations performed 
using formula (2) and the direct integration of the experimental data over the energy range studied 
in this work. As was explained above, the bare cross section was used in the calculation. The 

method is straightforward and has been described elsewhere [3]. The hrst line of the Table (Old data) 
gives the result based on the data of OLYA, CMD and DM1 while the second one (New data) is 
obtained from the CMD-2 data only. The third line (Old -|- New) presents the weighted average of 
these two estimates. The assumption about the complete independence of the old and new data used 
in the averaging procedure seems to be well justihed. For convenience, we list separately statistical 
and systematic uncertainties in the second column while the third one gives the total error obtained 
by adding them in quadrature. One can see that the estimate based on the CMD-2 data is in good 
agreement with that coming from the old data. It is worth noting that the statistical error of our 
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Table 5 


Contributions of the tttt channel to (g-2)/2 


Data 

a-V 10-10 

Total error, 10 1° 

Old 

374.8 ± 4.1 ± 8.5 

9.4 

New 

368.1 ± 2.6 ± 2.2 

3.4 

Old -|- New 

368.9 ± 2.2 ± 2.3 

3.2 


measurement is slightly larger than the systematic uncertainty. Because of the small systematic 
error of the new data, the uncertainty of the new result for is almost three times better than 
the previous one. As a result, the combined value based on both old and new data is completely 
dominated by the CMD-2 measurement. 

The example above only illustrates the importance of the improved accuracy. At the present time 
the analysis of the tttt data as well as other hadronic hnal states in the whole energy range accessible 
to CMD-2 is in progress. Independent information is also available [32] or expected in close future 
from other experiments studying low energy e+e“ annihilation [33]. When all the above mentioned 
data are taken into account, one can expect a signihcant improvement of the overall error of a!^'^ (by 
a factor of about 2) compared to the previous one based on the e+e“ data only [3]. 

The completion of the analysis of the tttt data will also provide a possibility of the precise test of the 
CVC based relation between the cross section of and the spectral function in the decay 

T~ — n~7i^h'r- The solution of the problem of the possible deviation between and r-lepton data 
[8,10] should also involve a thorough investigation of the effects of isospin breaking corrections as 
well as additional radiative corrections in r decays [21,34,35]. 


6 Conclusion 

The following values of the p and uj meson parameters have been obtained with the Gounaris-Sakurai 
ht to the formfactor data: 


M, 

= (776.09 ± 0.64 ± 0.50) MeV, 

rp 

= (144.46 ± 1.33 ± 0.80) MeV, 

r(p ^ e+e") 

= (6.86 ±0.11 ±0.05) keV, 

1 

+ 

T 

) = (1.33 ±0.24 ±0.05)%, 

arg(5 

= 12.6° ±3.7° ±0.2°. 


The measurement presented in this paper supersedes the preliminary result [16], obtained from the 
same data. Analysis of the much larger data sample collected by CMD-2 in 1996 (the 370-540 MeV 
energy range), 1997 (the 1040-1380 MeV energy range) and 1998 (the second scan of the 370-960 
MeV energy range) is in progress. 
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